Leaflet

Introduction to Leaflet in R

Basic Usage

Map Widget

create map widget with leaflet()

load in the data. We’ll use the beetle data again.

beetle_url <- "https://raw.githubusercontent.com/dyerlab/ENVS-Lectures/master/data/Araptus_Disperal_Bias.csv"
beetle <- read.csv(beetle_url)

beetle <- beetle %>%
   st_as_sf( coords=c("Longitude","Latitude"), crs=4326 )

Basemaps

Add a basemap with addTiles(). By default OpenStreetMap tiles are used.

Base maps are not required. Without a basemap your map will be floating in space. Using a basemap adds context, especially for point data

leaflet() %>%
  setView(lng = -77.5, lat = 37.5, zoom = 10) %>% #not necessary to set your view
  addTiles()

This is a blank map because we haven’t added any layers

Third party tiles

Use addProviderTiles() function to use a different basemap. Think about what information you are trying to convey when choosing a basemap. Will extra detail be helpful or will it distract from your data? Most third party base map styles have a labeled and an unlabeled option.

Base maps can be combined if, for example, you like the style of one basemap but the labels of another basemap. You can also include multiple basemap layers and let the viewer select their prefered map. We’ll see how to do this later.

See here for all options

leaflet() %>%
  setView(lng = -77.5, lat = 37.5, zoom = 10) %>% #not necessary to set your view
  addProviderTiles(providers$CartoDB.Positron)

Data Objects

Both leaflet() and each map layer have a data = parameter. Spatial data can be in the form of:

• dataframe with lat/lng columns
• sp objects
• sf objects
• maps package map() objects

Data can be passed through the leaflet() function or through the map layers.

Here we define our data in the leaflet() function

leaflet(beetle) %>%
  addTiles() %>%
  addCircles()

Here we define the data in the map layer (addCircles())

leaflet() %>%
  addTiles() %>%
  addCircles(data = beetle)

Point Data

Points can be plotted using addMarkers() or addCircles()

• Markers stay the same size regardless of zoom level

• Circles scale with the map

Markers

Default marker:

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addMarkers()

zoom in/out and notice that the markers stay the same size

Custom Markers

• Custom markers can be made using a url/image file or from libraries

• It is important to define the height and width (in pixels) of your icons

bug <- icons(
  iconUrl = "https://www.pngfind.com/pngs/m/14-144860_beetle-bug-png-transparent-image-bug-png-png.png",
  iconWidth = 20,
  iconHeight = 20
)

Use your custom marker in the icons = of the addMarkers() function

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addMarkers(icon = bug)

Marker from a library

link to Font Awesome library

makeAwesomeIcon() allows you to select your icon and edit its color ect.

addAwesomeMarkers() allows you to add your custom icon to the map

fa_bug <- makeAwesomeIcon(icon = "bug", library = "fa", 
                          markerColor = "cadetblue", iconColor = "beige")

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addAwesomeMarkers(icon = fa_bug)

NA

Notice that our custom icon using makeAwesomeIcon() and addAwsomeMarkers() displays the icon ontop of a marker. The custom icon using addMarkers() displays directly on the map

Popups and Labels

• Popups appear/disappear when clicked

• Labels appear/disappear when you mouse hovers on the object

Popups

Popups can be added as a stand-alone feature using addPopups(), or add to appear when a shape is clicked

Stand-alone

#stand-alone

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  setView(lng = -111, lat = 25, zoom = 7) %>%
  addPopups(lng = -111, lat = 25, paste0("beetles live here"))

NA
NA

Notice that the popup is its own leaflet layer

As a marker option

#as a part of a marker

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addAwesomeMarkers(icon = fa_bug, popup = paste0("Site:", beetle$Site))

NA

Notice that this time we define our popup inside the maker layer, not as its own layer

Labels

Use label = to add a label displayed on a mouse over

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addAwesomeMarkers(icon = fa_bug, label = paste0("Site:", beetle$Site))

Labels and Popups can be customized using the labelOptions = labelOptions() function. You can define the text size, color, font, box border, box shadow, ect.

addCircles()

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles()

Change radius:

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(radius = ~Males *500, stroke = FALSE, fillOpacity = .5)

NA

Change color:

pal <- colorNumeric(
  palette = "RdBu",
  domain = beetle$MFRatio,
  reverse = TRUE
)


leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(color = ~pal(MFRatio), fillOpacity = .7, radius = 15000, stroke = FALSE)

Highlight

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(fillColor = ~rev(pal(MFRatio)), fillOpacity = .7, radius = 15000, 
             weight = 1, color = "grey",
  highlightOptions = highlightOptions(color = "green", weight = 2,
      bringToFront = TRUE))

Map Groups

Assign layers to groups using group =. Assigning groups allows you to show/hide layers or control the visibility of layers through the function addLayersControl()

A group can be made up a a single layer or multiple layers, but each layer can only belong to one group.

leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(data = (beetle %>% filter(Males > Females)), 
             fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", radius = 10000, weight = 1,
             group = "males") %>%
  addCircles(data = (beetle %>% filter(Males < Females)), 
             fillColor = "pink", fillOpacity = 1, 
             color = "red", radius = 10000, weight = 1,
             group = "females") %>%
  addLayersControl(overlayGroups = c("males", "females"),
                   options = layersControlOptions(collapsed = FALSE))

NA

Using map groups with Provider Tiles

overlayGroups can be individually checked or unchecked

badeGroups can only be viewed one group at a time

leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron, group = "Default") %>%
  addProviderTiles(providers$OpenStreetMap, group = "Open Street Maps") %>%
  addProviderTiles(providers$Stamen.Toner, group = "Stamen Toner") %>%
  addCircles(data = (beetle %>% filter(Males > Females)), 
             fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", radius = 10000, weight = 1,
             group = "males") %>%
  addCircles(data = (beetle %>% filter(Males < Females)), 
             fillColor = "pink", fillOpacity = 1, 
             color = "red", radius = 10000, weight = 1,
             group = "females") %>%
  addLayersControl(overlayGroups = c("males", "females"),
                   options = layersControlOptions(collapsed = FALSE),
                   baseGroups = c("Default", "Open Street Maps", "Stamen Toner"))

Using map groups to control zoom levels with groupOptions()

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", radius = 10000, weight = 1,
             group = "beetles") %>%
  groupOptions("beetles", zoomLevels = 5:8)

Zoom in and out to watch the layer turn on/off

Clustering

When there are a larger number of makers on a map you can cluster them using clusterOpions =

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", weight = 1,
             group = "beetles",
             clusterOptions = markerClusterOptions()) 

Hovering over a cluster marker with your mouse allows you to see the coverage of the cluster

Click on a cluster to zoom in to the cluster bounds

All of these features can be changed markerClusterOptions. markerClusterOptions also allows you to freeze the clustering at a defined zoom level with freezeAtZoom =

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", weight = 1,
             group = "beetles",
             clusterOptions = markerClusterOptions(
               showCoverageOnHover = FALSE,
               freezeAtZoom = 6
             )) 

Now the clusters stay the same regardless of zoom level. When the cluster is clicked you can see each point included in the cluster

Legends

Legends are added using addLegend()

• must provide color information through pal = if using a palette, or color = if using custom defined colors

• must provide the values used to generate colors from the palette or the labels for each corresponding color

leaflet() %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(data = (beetle %>% filter(Males > Females)), 
             fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", radius = 10000, weight = 1,
             group = "males") %>%
  addCircles(data = (beetle %>% filter(Males < Females)), 
             fillColor = "pink", fillOpacity = 1, 
             color = "red", radius = 10000, weight = 1,
             group = "females") %>%
  addLegend( colors = c("lightblue", "pink"),
             labels = c("Males", "Females"),
             opacity = 1)

NA

Example of legend with continuous data

pal <- colorNumeric(
  palette = "RdBu",
  domain = beetle$MFRatio,
  reverse = TRUE)

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircles(fillColor = ~pal(MFRatio), fillOpacity = .7, radius = 15000, 
             weight = 1, color = "grey") %>%
  addLegend(pal = pal, values = ~MFRatio,
            title = "Male to Female Ratio")

Scale Bars

Scale bars are added using addScaleBar()

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", weight = 1) %>%
  addScaleBar(position = "bottomright")

watch the scale bar adjust itself as you zoom in and out

Leaflet Extras Package

The leaflet.extras package allows for additional functionality via leaflet plugins. We’ll take a look at three examples here, but there are tons of plugins to choose from

addResetMapButton() resets you map to the original view and zoom level

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", weight = 1) %>%
  addResetMapButton()

Move the map and zoom in/out, then press the reset button to watch the map return to its starting view and zoom level

addSearchOSM() allows you to search for locations on the map

library(leaflet.extras)

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", weight = 1) %>%
  addSearchOSM(options = searchOptions(
    zoom = 8,
    hideMarkerOnCollapse = TRUE,
    autoCollapse = FALSE
  ))

Try searching for “Cabo San Lucas”

addDrawToolbar() allows you to draw points, lines and polygons on map

leaflet(beetle) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(fillColor = "lightblue", fillOpacity = 1, 
             color = "blue", weight = 1) %>%
  addDrawToolbar(
    editOptions = editToolbarOptions(
      selectedPathOptions = selectedPathOptions()
    )
  )  %>%
  removeDrawToolbar()

NA

Try drawing and then deleting a shape from your map

Making a Choropleth Map

Lines and Polygons work in a very similar manner to circles

For this example we’ll use census tracts downloaded using the tigris package

library(tidyverse)
library(tigris)
options(tigris_use_cache = TRUE)

tracts <- tracts("VA", "Richmond city") %>%
  st_transform(4326)

Using FIPS code '51' for state 'VA'
Using FIPS code '760' for 'Richmond city'

Defining the color palette using the land area of each tract. Area is given in square meters. To convert to square miles we divide our area by 2.59e+6

pal <- colorNumeric(palette = "Blues",
                    domain = tracts$ALAND/2.59e+6) 

The polygon colors, which define a choropleth map, are assigned by fillColor = ~pal(ALAND/2.59e+6)

Labels, a legend and a scale bar are added to give the map a more polished finish

leaflet(tracts) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addPolygons(fillColor = ~pal(ALAND/2.59e+6),  fillOpacity = 1, smoothFactor = .2,
              color = "darkgrey", weight = 1,
                highlightOptions = highlightOptions(color = "white", weight = 2, opacity = 1,
                bringToFront = TRUE),
              label = paste(round(tracts$ALAND/2.59e+6, digits = 2), "sq. mi.")
              ) %>%
  addLegend(pal = pal, values = ~ALAND/2.59e+6, title = "Area (sq. miles)") %>%
  addScaleBar(position = "bottomright")

NA

The map is now shaded so that the larger census tracts are darker blue and the smallest tracts are lighter blue